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Abstract 

Optical crosstalk seriously limits the photon-counting resolution of silicon photo- 
multipliers. In this work, realistic analytical models to describe the crosstalk effects 
on the response of these photodetectors are presented and compared with experi- 
mental data. The proposed models are based on the hypothesis that each pixel of 
the array has a finite number of available neighboring pixels to excite via crosstalk. 
Dead-time effects and geometrical aspects of the propagation of crosstalk between 
neighbors are taken into account in the models for different neighborhood configura- 
tions. Simple expressions to account for crosstalk effects on the pulse-height spectrum 
as well as to evaluate the excess noise factor due to crosstalk are also given. Ded- 
icated measurements were carried out under both dark-count conditions and pulsed 
illumination. Moreover, the influence of afterpulsing on the measured pulse-height 
spectrum was studied, and a measurement of the recovery time of pixels was reported. 
High-resolution pulse-height spectra were obtained by means of a detailed waveform 
analysis, and the results have been used to validate our crosstalk models. 



1 Introduction 



Silicon photomultipliers (SiPMs) belong to the sort of photodetectors with single-photon 
detection capability [1J . These devices consist in a monolithic array of multiple single ava- 
lanche photodiodes (APDs) operated in Geiger mode. Each component or pixel produces 
a pulse of constant amplitude regardless of the number of impinging photons. All pixels 
are connected to the same output channel providing a total output signal equal to the 
sum of those from the individual pixels, which enables a large dynamic range and photon 
counting. SiPMs offer desirable properties like high gain and sensitivity, and very good 
time and photon-counting resolutions, which have increased interest in the use of these 
photodetectors in numerous applications displacing more traditional ones. Unfortunately, 
their photon-counting ability is seriously limited by optical crosstalk, which affects the 
linearity of the detector, in particular, causing a significant excess noise. The very rele- 
vant influence of crosstalk on the performance of SiPMs has been contemplated in most 
characterization studies of these devices. However, a more accurate description of this 
effect is still necessary. 

When a pixel is fired by either an incoming photon or by a thermally generated electron- 
hole pair, hot carriers in the avalanche breakdown induce emission of IR photons [2HI] that 
in turn may trigger further avalanches in nearby pixels. This stochastic process, called 
optical crosstalk, is characterized by being nearly instantaneous, and its probability is 
proportional to the SiPM gain. Operating at low bias voltage would diminish significatively 
crosstalk effects, but at the expense of degrading the photon-detection efficiency. The 
incorporation of isolation trenches around each pixel, as proposed by [5], successfully 
reduces optical crosstalk [6HH]- This technique has become very usual in the fabrication 
of this kind of photodetectors. However, it was first shown in [5] that a very significant 
contribution to crosstalk can come from light reflected on the bottom surface of the Si 
bulk, and thus, trenches would not be able to completely prevent it. 

The main effect of crosstalk is to introduce a multiplication noise, although it does 
not affect the pulse-height resolution. So, at conditions where only one pixel is expected 
to be excited simultaneously (e.g., dark counts), crosstalk results in output pulses with 
amplitudes twice or several times the amplitude of a single triggered pixel. The crosstalk 
probability e is usually defined as the rate of dark counts with crosstalk (two or more fired 
pixels) divided by the total dark-count rate. 

Several studies based on Monte Carlo simulations [101412] have shown that the response 
of various SiPM devices from Photonique and Hamamatsu can be properly described by 
assuming that crosstalk only takes place between adjacent pixels. On the other hand, 
complete simulations [13] of some non-commercial SiPM devices, including light propa- 
gation through the silicon bulk and reflections on the bottom surface, have demonstrated 
that crosstalk can also be induced in distant pixels, which is supported by experimental 
data for these SiPMs. 

An analytical formulation is desirable for a better understanding of crosstalk effects. 
Most of the available statistical models of crosstalk [14H21] are based on simple assump- 
tions involving restrictions on the number of crosstalk events per initially fired pixel or on 
crosstalk cascading (i.e., crosstalk events generated by other crosstalk-induced avalanches). 
S. Vinogradov has developed very recently an analytical model overcoming these limita- 
tions, where crosstalk is considered as a branching Poisson process [22]. However, no 
analytical model has so far dealt with the geometrical arrangement of the neighboring 
pixels available to be excited via crosstalk and the consequent local saturation effects as 
pixels become inactive after being excited. In the present work, realistic analytical mod- 
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els of crosstalk that take into account these geometrical considerations are reported and 
validated by comparison with experimental data. 

This paper is organized as follows. In section [2j the proposed crosstalk models are 
explained, and analytical expressions are provided for the probability distribution of the 
number of crosstalk events at both dark-count and pulsed-illumination conditions. The ex- 
perimental setup and the analysis method are described in section [3j Results are presented 
in section |U Conclusions are drawn in section 03 

2 Analytical models of crosstalk 

Next, fundamentals of the statistics of crosstalk events (including cascading effects) for 
a single initially fired pixel are described. Then, realistic crosstalk models that take 
into account the above-mentioned geometry-dependent saturation effects are detailed in 
section [2~2) Finally, effects of crosstalk on the photon statistics under pulsed-illumination 
conditions are discussed in section l2~3l 

2.1 Fundamentals 

Several analytical models of crosstalk noise in SiPMs are available in the literature. Many 
of them [14H19j assume that crosstalk obeys a Bernoulli distribution, that is, a primary 
avalanche can either trigger a secondary avalanche in a neighboring pixel with probability 
p or no avalanche with probability 1 — p. Some of these models also include cascading 
processes where a secondary avalanche may trigger a tertiary avalanche, which could in 
turn trigger a quaternary avalanche and so on, although the number of crosstalk events 
that can be induced directly by each pixel is limited to one. However, a single avalanche is 
made up of a large number of carriers (typical avalanche gain is 10 5 - 10 6 ), each one being 
able to induce emission of crosstalk photons with a probability of ~ 3 • 10 -5 [31 [131 122]. 
Accordingly, some authors [20-22J have used a Poisson distribution instead to describe the 
number of crosstalk events initiated directly by a single avalanche. 

Even though many IR photons may be produced in a single avalanche breakdown, they 
are expected to be absorbed preferably in the vicinity of the pixel in which the avalanche 
was triggered, and therefore, the number of neighboring pixels able to be excited efficiently 
via crosstalk must likely be small. This has suggested to us to use a binomial distribution 
to describe the number s of successful crosstalk events in a finite sample of n available 
neighbors of the primary pixel, that is, 

P{s)= Qp s (l- P ) n - S , (1) 

where the simple approximation of assuming the same probability p of crosstalk for any 
individual neighbor is made. This p parameter is related to the total crosstalk probability 
e (one or more crosstalk events) through 1 — e = (1 — p) n . Note that ([I]) reduces to 
the Bernoulli distribution for n = 1 and tends to a Poisson distribution with mean A = 
— ln(l — e) for n — > oo keeping constant e. The latter case would correspond to the limit 
situation where crosstalk photons originated from the primary pixel can reach any other 
pixel of a very large array. 

To account for the total number of crosstalk events per initially fired pixel, crosstalk 
cascading should be included. As a first approximation, we assume that any triggered pixel, 
either the primary pixel or a secondary or higher-order one, can always induce crosstalk 
in n more neighboring pixels. So, the probability of a cascade of crosstalk events would 
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Figure 1: Histories of crosstalk excitations for a model of 2 neighbors and total number of 
triggered pixels from 1 to 3. The primary and the subsequent crosstalk events are represented by 
grey-filled and red-filled circles, respectively. The cascade development is indicated by solid arrows, 
while a dotted arrow ending in an open circle means that the cascade no longer continues because 
the neighbor is not triggered. The corresponding probabilities are given in terms of combinatorial 
numbers. 



be obtained by applying ([I]) repeatedly for each pixel of the sequence. The probability to 
get a certain total number k of triggered pixels is proportional to the number of histories 
with k — 1 crosstalk events, considering all the possible combinations of neighbors that are 
excited or not in each step. This is illustrated in figure [TJ where all the histories with up 
to two crosstalk events are identified for a model of n = 2 neighbors. For arbitrary n and 
k values, probabilities can be calculated as 

P(k) = h n ^ 1 p k - 1 (l-p) kn - k+1 , (2) 

where /i n ,fc-i denotes the number of histories of crosstalk excitations adding up to k — 1 
for a model of n neighbors. The count of histories is not trivial for a large number of 
triggered pixels. Nevertheless, the procedure shown in figure [T] can be generalized by the 
recursive formula 
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hn,c = ^ ] ( j ^ ] h n i L ^ ^ h n ^ 2 . . . ^ ] h n i s _ 1 h n)C - s -i 1 -.„-i s _ 1 , (3) 

s=l ^ ' j 1= i 2 =0 « s -i=0 

starting with h n Q = 1. This results in the following probabilities for k < 5: 



P(l) = (l-p)" = l- £ 
P(2) =np(l-p) 2n - 1 

P(3) = ^n(3n-l)p 2 (l-p) 3n ~ 2 

P(4) = l n (8n 2 - 6n + l) p 3 (l - p) 4n ' 3 
3 

P( 5 ) = \J l^ n 3 _ 25 n 2 + —n - l) p\l - pf n ~ A . (4) 
4 \ 6 6 J 

Notice that the overall crosstalk probability e is unaffected by including or not cascading 
effects, because it only relies on the original binomial probability that the primary pixel 
induces no crosstalk event. 
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Figure 2: Probability distributions of the total number of triggered pixels for a single initially 
fired pixel resulting from equations ^ and (|3| for different numbers of neighbors n. The total 
crosstalk probability (one or more crosstalk events) is set to e — 0.5 (left) and to e = 0.15 (right). 



The mean and variance of the distribution are conveniently approximated as power 
series of either p or e up to second ordeiQ 

E = 1 + np + (np) 2 + o (p 2 ) = 1 + e + 3n ~ 1 e 2 + o (e 2 ) (5) 

' In 

7n — 3 

Var = np + n(3n — l)p 2 + o (p 2 ) = e H e 2 + o (e 2 ) . (6) 

2n 

In figure [2j results of equations ^ and ^ are shown for two different crosstalk prob- 
abilities e = 0.5 (left) and e = 0.15 (right) and for several n values. The probability 
distributions for n = 1 and for n — > 00 correspond respectively to the geometric and Borel 
distributions, which have already been exploited by S. Vinogradov et al. to model crosstalk 
|16| I22|. For intermediate n values, the probability distribution is between these two limit 
situations, although it tends rapidly to the Borel distribution as n increases. In general, 
the larger the number of neighbors is, the more likely multiple crosstalk excitations are. 

Crosstalk can be considered as an extra multiplication process introducing noise. This 
is usually characterized by the excess noise factor according to the standard definition 

\f (IT 

ENF = 1 + ^, (7) 

where E represents the mean multiplication factor (gain) of the signal due to crosstalk for 
a single initially fired pixel. Using ([5]j6]), the excess noise factor for a crosstalk model of n 
neighbors can be approximately expressed as 

3n — 3 

ENF = l + np + n(n - l)p 2 + o (p 2 ) = 1 + e H e 2 + o (e 2 ) . (8) 

2n 

This parameter increases as n increases at constant e, as expected, although it saturates 
at a value of 1/[1 + ln{\ — e)], which corresponds to the Borel distribution. 

2.2 Models including saturation effects 

The frame of a binomial crosstalk probability with n neighbors per pixel is our starting 
point to construct a realistic crosstalk model. Probability distributions resulting from 



1 Both mean and variance grow rapidly, and even diverge, at increasing p and e values. As a consequence, 
this approximation fails at large crosstalk probability (e > 0.2). 
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Figure 3: Neighborhoods defined in the proposed geometrical crosstalk models including satu- 
ration effects. From left to right: "4 nearest neighbors", "8 nearest neighbors", "8 L-connected 
neighbors" and "all neighbors" . 



account for cascading of crosstalk; on the other hand they do not include 
saturation effects due to the fact that a pixel that has already been triggered is inactive 
during a short period of time. 

Different pixels can share neighbors in such a way that the last pixels in a history of 
excitations have less chances to induce further crosstalk events. In fact, if pixel i is a 
neighbor of pixel j, pixel j is also a neighbor of pixel i, but crosstalk between both pixels 
can only take place in one direction at a time: i — )■ j or j — >■ i. Thus, in a model of n 
neighbors, only the primary pixel actually has the full number of neighbors available to 
excite. This saturation effect will depend on how dense the neighborhood of a pixel is. To 
evaluate it in an analytical way, we have considered four geometrical models of crosstalk 
defining different neighborhoods in a square lattice of pixels: "4 nearest neighbors", "8 
nearest neighbors", "8 L-connected neighbors" and "all neighbors" (figure [3]). 

The 4-nearest-neighbors hypothesis has already been used successfully in Monte Carlo 
simulations of the response of SiPMs including the crosstalk effect |1UH12| . However, as 
pointed out above, some works have shown that crosstalk is also possible, or even more 
probable, between non-contiguous pixels [9l [13], which has motivated to evaluate the 
other models. Note that the range of IR photons reflected on the surface of the detector 
is expected to be strongly dependent on the particular chip configuration (e.g., pixel 
separation, width of the Si bulk and depth of isolation trenches), and therefore, crosstalk 
effects in different SiPMs can be better described by different geometrical models. 

Unfortunately, simple expressions of the probability distribution of arbitrary number 
k of triggered pixels cannot be derived for such models. Nevertheless, only a few crosstalk 
excitations per initially fired pixel are likely to occur at moderate crosstalk probability, 
and it may suffice, in practice, to determine the first four or five probabilities of the 
distribution. The probability of k = 1 (i.e., no crosstalk event) is trivially given by 
P(l) = (1 — p) n = 1 — e. To calculate the probabilities of k > 2, one can proceed as 
follows: 

1. All the different histories with k— 1 crosstalk events are identified for the actual model 
taking into account the geometrical arrangement of neighbors and pixel inactivation. 

2. The probability of each history is calculated as p c (1 — p)f , where c = k — 1 is the 
number of crosstalk events and / is the overall number of crosstalk fails considering 
the active neighbors left for each triggered pixel following excitation order. 

3. The final probability P{k) is computed as the sum of probabilities of all the con- 
tributing histories. 



equations (2p 
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Figure 4: Several crosstalk histories (k = 5) giving rise to the same pattern of excitations in the 
4-nearest-neighbors model. Both the number c of crosstalk excitations (arrows) and the number / 
of crosstalk fails (short lines ending with dots) are indicated. Excitation order (see labels on the 
corner of triggered pixels) is followed to determine the active neighbors left for each pixel. To avoid 
ambiguities, crosstalk excitations induced by the same parent are assumed to be ordered clockwise 
starting from the top pixel. 

Table 1: Probabilities of k < 5 triggered pixels for the different proposed models. The q parameter 
is defined as 1 — p and N is the total number of pixels in the array (last column) . 
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For instance, assuming the 4-nearest-neighbors model, there are four symmetrical his- 
tories with k = 2, where the primary pixel and one of its neighbors are triggered, but nei- 
ther the remaining three neighbors of the primary pixel nor the other three neighbors left 
for the secondary one are excited by crosstalk (i.e., / = 6). Therefore, P(2) = 4p(l — p) 6 
is obtained. For k = 3, 18 histories can still easily be identified, all of them with / = 8, 
resulting in P(3) = 18p 2 (1 — p) s . Calculations become involved for larger k values, since 
the number of histories increases rapidly, and they can also have different number / of 
crosstalk fails due to the fact that pixels share neighbors. To illustrate this, a few histories 
with k = 5 are shown in figure |4j 

A computer algorithm was made to sequence the histories with k < 5 for the 4-nearest- 
neighbors, 8-nearest-neighbors and 8-L-connected models. An endless array of pixels was 
assumed, that is, border effects were ignored. Collecting histories with same k value, 
analytical expressions of the corresponding P{k) probabilities were obtained for these 
models. For the all-neighbors model, probabilities of k < 5 were also determined as a 
function of the p parameter and the number iV of pixels in the array using combinatorial 
mathematics. Results are listed in table [TJ 

Figure [5] displays the probability distributions up to k = 5 calculated for the 4-nearest- 
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Figure 5: Probability distributions for the 4-nearest-neighbors model (left) and for both the 8- 
nearest-neighbors and 8-L-connected-neighbors models (right) with e = 0.4. Results are compared 
with the distributions of equations 
saturation effects. 



2p| for n — 4 and n = 8, respectively, to illustrate the 



neighbors model (left) and for both the 8-nearest-neighbors and 8-L-connected-neighbors 
models (right) setting e = 0.4. The comparison with the distributions of equations ^ 
[3]) for n = 4 and n = 8, respectively, shows a redistribution of probabilities due to 
saturation, probabilities of k = 2 and k = 3 being increased while those of more triggered 
pixels decreased. This can be interpreted as a consequence of the reduction of the effective 
number of active neighbors (see figure]!]). Saturation effects are significant in the 4-nearest- 
neighbors model, while they are rather small and weakly dependent on the geometrical 
arrangement of the neighborhood in both models of 8 neighbors. In fact, saturation is 
almost negligible in the all-neighbors model for N > 100, which yields results very similar 
to the Borel distribution. 

The mean, variance and excess noise factor of the probability distributions given in 
table [1] were found to be described, up to second order of approximation in p or s, by the 
following unified expressions: 

3n — 3 

E = l + np + n(n- l)p 2 + o (p 2 ) = 1 + e H e 2 + o (e 2 ) (9) 

v ' 2n 

7n — 9 

Var = np + n(3n - 4)p 2 + o (p 2 ) = e H e 2 + o (e 2 ) (10) 

2n 

3n — 9 

ENF = l + np + n(n- 4)p 2 + o (p 2 ) = 1 + e H e 2 + o (e 2 ) , (11) 

In 

only depending on the number of neighbors assumed in each model, i.e., n = 4, 8 or N — 1. 
In a third or higher-order approximation, this is no longer possible since the particular 
geometrical arrangement of the neighborhood also plays a role in the development of 
crosstalk histories. It should be noted that the E, Var and ENF parameters for the 
proposed geometrical models are lower than those for the models of section |2.1| for the 
same number of neighbors but without saturation effects (see equations (J5]|8j) ) . 

We checked that a geometric extrapolation is a good approximation for the probabilities 
of k > 5 for the proposed geometrical models 



P(k) ~ P(5) 
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where the geometric rate of decrease is set so that X^bLi -P(&) = 1- For instance, this 
allows us to compute numerically the mean, variance and excess noise factor of the above 
probability distributions for larger e values and with better precision than equations ([911). 



2.3 Photon-counting statistics under pulsed illumination 

The above probability distributions are derived assuming a single initially fired pixel and 
they can be applied to describe the response of SiPMs with crosstalk at dark conditions 
or low-intensity illumination. For measurements under pulsed illumination, simultaneous 
random excitation of primary pixels (usually Poisson distributed) has to be accounted for. 
The composition of a Poisson distribution of initially fired pixels with either a geometric 
or a Borel distribution (cases n = 1 and n — > oo in equations ([2]{3])) has already been used 
by \16\ I22j. In a general case, the probability distribution of the total number of triggered 
pixels, k = m primaries +c crosstalk events, is given by 



Ptot(O) = Ppr(O) 
k 

Ptot{k) = PpAm) P m {k) , k = l, 2, 3,..., (13) 

m=l 

where P pr (m) is the probability distribution of the number m of primary events and P m (k) 
is the probability of total number k of triggered pixels provided m primaries. The latter 
probability can be calculated from the recursive formula 

k—m 

P m +i(k) = Pm(k-i)Pi(i), m = l, 2, 3,..., (14) 

i=l 

with Pi(i) being the probability distribution of the number i of triggered pixels for a single 
initially fired pixel, like those given in table [T] for the proposed geometrical models. 

These expressions do not account for possible further saturation effects caused by 
primary events in close pixels initiating crosstalk histories that overlap each other. Nev- 
ertheless, if the number of primary events is low compared to the total number of pixels 
of the device, the probability of overlapping crosstalk histories is negligible and this ap- 
proximation can be done. In fact, more than 20% of pixels of the array should be fired 
initially so that this saturation effect is important for typical crosstalk probabilities in the 
4-nearest-neighbors model |12| . 

In principle, the original photon statistics P pr (m) can be extracted from the observed 
probability distribution Ptot(^) if the crosstalk distribution P\{k) is known (e.g., from the 



pulse-height spectrum of dark counts). Using (13 14), the following sequence is derived: 



Ppr(0) = Ptot(0) 

p m - Ptot ^ 

P pT (m) 



(l-e)< 



m—l 



PtotM- P pi (l) Pi(m) 



m = 2, 3, 4, 



(15) 



l=i 
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However, this can more easily be accomplished making simple assumptions on the expected 
photon statistics. In the most usual case of a Poisson distribution, the photon statistics 
is entirely determined by the mean, which is simply given by 



H= -ln(P tot (0)) , (16) 

where -Ptot(O) is unaffected by crosstalk indeed. Once \i is obtained, the crosstalk proba- 
bility can be calculated as 

s = l- Pto f\ . (17) 
ii exp(-/i) 

Note that expressions ( 16 ) and ( 17 ) are independent of the crosstalk model, and we will 
use them in section [4] to calculate respectively the fi and e parameters from experimental 
dataj^J In other cases, however, it may not be possible to use this approach either because 
the photon arrival is not Poisson distributed or because the -Ptot(O) and Ptot(l) probabilities 
cannot be determined accurately (e.g., if the number of impinging photons is large, these 
probabilities will be too small). 



From (13 14), it can be demonstrated that the mean and variance of the observed 



probability distribution are related to those of photons and crosstalk events as 

^tot = E pT Ei (18) 

Vartot = E pv Van + Var pr E\ , (19) 

where E\ and Var\ can be calculated using the approximate expressions given in table [T] for 
the proposed crosstalk models)^] Therefore, valuable information of the original probability 
distribution of photons can be obtained simply from the mean and variance of the measured 
pulse-height spectrum if the crosstalk probability is known. 

To evaluate the excess noise introduced by crosstalk, the contribution of the random 
excitation of primary pixels (i.e., input noise) should be discounted from the total noise 
in the output signal. Following |16l \2'1\ , the excess noise factor due to crosstalk can be 
defined as the relative losses in the signal-to-noise ratio from input to output 

ENF> = E } IVa ^ . (20) 
E 2 tJVart t K ' 

However, this definition is not equivalent to the standard one given by (|7|) for dark-count 
measurements. We found that the standard definition can be generalized for pulsed- 
illumination conditions using the following alternative expression: 

Vartot Var pr \ 1 , Var\ 



E ^ l+ M^r-^r^ l+ ^p <2l) 

Note that when E pr = Var pr , as happens for a Poisson distribution, both expressions 
(20) and (21) are identical, whereas only (21) properly reduces to Q when E pi = 1 and 



Var pr = (e.g., dark counts with a single primary event). 



2 This procedure has previously been employed in [llj . 

3 In this paper, we use a simplified notation in which the names of the random variables are omitted 
(symbols for their numerical values are only shown), and indexes are used instead to indicate the cor- 
responding random variable. For instance, Ptot(fc), -Etot and Var to t stand respectively for P(A"tot = k), 
E[K to t] and l/ar[_K'tot], where Ktot is the random variable associated to these probability distribution, 
mean and variance. 
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Figure 6: Schematic view of the experimental setup consisting of a SiPM (Hamamatsu S10362- 
11-100C) with associated electronics, a N2 laser emitting 3.5 ns FWHM pulses and a digital oscil- 
loscope. 



3 Experimental method 

A dedicated experiment was carried out to study the crosstalk effect in a SiPM device. 
Two types of measurements were performed: under low-light-level conditions (dark counts) 
and using a pulsed light source to illuminate the device. 

The experimental setup is described in section [3~T| The data analysis, which includes 



the processing of signals with a specific software, is explained in section 3.2 



3.1 Setup 

A schematic view of the experimental setup is shown in figure [6j This includes the SiPM 
with the associated electronics, the light source and the data-acquisition system. The 
SiPM chosen for this study was the model S10362-11-100C from Hamamatsu with 100 
pixels in an active area of 1 mm x 1 mm. It was connected to a bias circuit that had 
a load resistor of 1 kO and a coupling capacitor of 0.1 fj,F. Signals were sent to a fast 
amplifier based on AD8367 (Analog Devices) with a bandwidth of 500 MHz and nominal 
gain of 42.5 dB. The amplifier was configured to provide narrow output pulses with a rise 
time of about 2 ns and an exponential decay with 23 ns time constant. Typical pulse sizes 
of 20 mV were achieved for single-photoelectron signals when working at an overvoltage 
of around 1.2 V. 

For the measurements under pulsed illumination, we used a 337 nm nitrogen laser 
(Stanford Research Systems, model NL100) emitting light pulses of 3.5 ns FWHM at a 
repetition rate of up to 20 Hz. The laser beam was attenuated and arranged to illuminate 
uniformly the active area of the detector. The laser-pulse intensity was monitored using 
an internally installed photodiode, which also provided the synchronized output signal. 

Data were recorded using a digital oscilloscope (Tektronix TDS5032B) with a band- 
width of 350 MHz and a sampling rate of up to 5 GS/s. The oscilloscope enabled us 
to store all the time-intensity information of signals in ASCII files for later analysis. In 
dark-count measurements, each file contained a 400 fj,s long run, while in pulsed regime, 
the data acquisition was done within a time window of 200 ns in coincidence with the laser 
pulse (4000 laser shots per file). The time resolution was set to 0.4 ns per channel in both 
cases. 

Both bias voltage and temperature were monitored. The bias voltage was stable within 
0.1% during each run. Temperature was measured using a type K thermocouple close to 
the SiPM, enabling us to register changes of 0.5 °C. 
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Figure 7: Dark counts registered at an overvoltage of 1.2 V. Left: Identification of pulses through 
the signal deconvolution. Right: Determination of the pulse amplitude using an exponential fit to 
determine the baseline of each pulse. 



3.2 Waveform analysis 

A C++ program was developed to analyze the raw data stored by the digital oscilloscope. 
Essentially, the algorithm searched for pulses and then measured their arrival times and 
amplitudes. A linear digital filtering of the signal was also performed to reduce the elec- 
tronic noise. Some sample pulses registered at dark conditions are shown in figure [7] (the 
digital filter is already applied to signals). 

Pulses were identified by their sharp leading edge. To this end, the signal was previ- 
ously deconvolved with an exponential decay function with 23 ns time constant, in such 
a way that tails of pulses were suppressed and edges were clearly distinguished even for 
piled- up pulses (left-hand plot of figure [7]) . Threshold criteria were applied to this decon- 
volved signal on both amplitude and pulse width to further discriminate electronic noise. 
Then, the arrival time of an identified pulse was measured as the time at which the de- 
convolved signal reaches the maximum, corresponding to the position of the pulse edge 
in the original signal. The smallest time difference to resolve close pulses was determined 
to be 3.0 ± 0.6 ns, which was our coincidence time window in dark-count measurements 
(crosstalk is assumed to be nearly instantaneous). On the other hand, when illuminating 
with the laser beam, light was detected within a time interval of 5 ns, which was selected 
as the coincidence window instead, that is, two pulses that arrive within this time interval 
are considered as a single pulse even though they can actually be resolved. 

The pulse amplitude was measured as the peak height from the baseline. The original 
signal was used for this purpose, because using the deconvolved signal resulted in larger 
uncertainties. Since pulses are often on the tail of the preceding one (see figure [7]), the 
baseline of each individual pulse was calculated by fitting the function 

f(t)=A + Bexp(-~)> with r = 23 ns, (22) 

to data in the time interval of 8 ns prior to the pulse, and extrapolating it up to the 
peak time. This method was proved to significantly enhance the measured pulse-height 
resolution of both dark counts and laser pulses. For instance, the amplitudes of the pulses 
of the right-hand plot of figure [7] were determined to be 44.5 mV and 20.3 mV, which 
are typical amplitudes of pulses having two and one triggered pixels, respectively, at an 
overvoltage of 1.2 V. However, limitations were found for piled-up pulses. If two pulses 
are closer than 20 ns, the amplitude of the second one cannot be accurately calculated 
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because the fit of equation (22 ) to data fails. In addition, if the time difference is less than 
10 ns, the amplitude of the first pulse cannot be measured either, because the peak is not 
well distinguished from the second pulse. Therefore, although the piled-up pulses of the 
left-hand plot of figure [7] are resolved (time difference of 6.4 ns), their amplitudes would 
not be measured separately. Quality cuts were applied to the identified pulses to avoid 
these problems, resulting in a further improvement of the pulse-height resolution. For the 
evaluation of crosstalk, a dead time of 120 ns was also introduced after each identified 
pulse to prevent afterpulsing effects (see section 4.1). 

Our analysis software generated a table of arrival times and amplitudes of the identified 
pulses classified according to the quality cuts and dead-time criteria. This allowed us to 
obtain pulse-height spectra and to study possible time-amplitude correlations as well as 
the effect of the applied cuts. The algorithm also calculated the rate of pulses in dark- 
count measurements. For this purpose, the acquisition time was calculated taking into 
account the applied dead time and quality cuts. 



4 Results 



Experimental results are presented and discussed below. In the first place, amplitude 
and timing measurements of afterpulsing, relevant for a correct treatment of experimental 
data, are reported. In section 4.2, the probability distribution of the number of triggered 
pixels is extracted from the pulse-height spectrum measured at dark-count conditions and 
it is used to discriminate between the crosstalk models described above. In section 14.31 



applicability of equations (18) and (19) to describe the effects of crosstalk on some pulse- 
height spectra measured under pulsed-illumination conditions is evaluated. 



4.1 Afterpulsing and recovery time 

Parasitic avalanches can be produced by the delayed release of carriers trapped in deep- 
level defects during a previous avalanche [24, 25J. This effect is referred to as afterpulsing 
and, unlike crosstalk, it occurs in the same pixel where the primary avalanche was devel- 
oped. After breakdown, a pixel need some time to recover the original bias voltage, so 
secondary avalanches with a delay less than the recovery time are weaker than ordinary 
ones. Short-delay afterpulses usually appear as small humps on the tail of the primary 
pulse and do not affect the pulse-height measurement. However, they can be counted as 
separate pulses of lower amplitude, distorting the pulse-height spectrum. Furthermore, 
afterpulsing avalanches are obviously able to induce crosstalk, but low-intensity ones are 
expected to have a smaller probability, and thus, they would contribute to reduce the 
average crosstalk probability if they are not discriminated. As pointed out in the previous 
section, our analysis algorithm applies a dead time long enough to prevent these effects 
due to afterpulse contamination. 

The recovery time of our SiPM was estimated using a technique similar to [26\ 127] . 
that is, studying the dependence of the afterpulse amplitude on the delay relative to the 
primary pulse. To this end, dark counts were registered continuously in a wide time win- 
dow and analyzed without imposing dead time so that afterpulses were also included. 
Given any two consecutive pulses, they may be either uncorrelated events or an after- 
pulse following its primary event. Uncorrelated avalanches are generally induced in fully 
recharged pixels, and thus, they should have constant amplitude, while afterpulses have 
an amplitude proportional to the overvoltage reached in the pixel since the last avalanche 
breakdown, with an expected exponential recovery. Both contributions are clearly visible 
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At (ns) 



Figure 8: Two-dimensional histogram of the pulse amplitude versus the arrival time relative to 
the preceding pulse. Two branches can be seen due to contributions of non-correlated events and 
afterpulses, respectively. An exponential fit to the lower branch is shown, providing an estimate 
of the recovery time. Crosstalk effects result in the repetition of the same pattern of branches at 
higher amplitudes. 



in the two-dimensional histogram displayed in figure [8j showing that pulse amplitudes 
are distributed in a narrow band that splits into two branches when the time difference 
with respect to the preceding pulse is smaller than ~ 80 ns, the lower branch being due 
to short-delay afterpulses. To better resolve these contributions, pulses closer than 80 ns 
to their two preceding pulses were rejected, that is, potential short-delay afterpulses that 
are not produced by the immediately preceding avalanche (fired in a different pixel) were 
eliminated. Crosstalk effects result in the repetition of the same pattern with two branches 
at higher amplitudes. The cut at 20 ns is due to the limitations of our analysis algorithm 



to calculate the amplitude of piled-up pulses (section 3.2). 

As expected, the afterpulse amplitude grows exponentially at increasing delay, although 
the recovery curve was found to be more suitably described if an offset time of 5 ns is added. 
We measured a recovery time (i.e., the reciprocal of the exponential rate) of 40 ns with 
an estimated uncertainty of 5 ns, which is consistent with the result reported in |26j for 
the same SiPM device and at room temperature. 

4.2 Dark-count measurements 

The pulse-height spectrum of dark counts was measured at an overvoltage of 1.2 V and 
room temperature (25 °C). In figure the experimental spectrum is shown with and 



without applying the cuts described in section 3.2 The controlled dead time introduced in 



the data analysis reduces significantly the statistics, while it also removes the pronounced 



shoulder on the left of each peak due to low-amplitude afterpulses (see section 4.1 ). Quality 
cuts discriminate a very small fraction of piled-up pulses that mostly populate the valleys, 
resulting in a slight improvement of the height-pulse resolution with no significant bias in 
the experimental probability distribution. Five Gaussian-like peaks were clearly resolved 
in the spectrum corresponding to pulses from 1 to 5 triggered pixels, although pulses of 
higher amplitudes were also registered. The number of collected events passing the cuts 
was 3.1 • 10 5 . 

The probability that a pulse has k = 1, 2, 3, 4, 5 triggered pixels, P(k), was determined 
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Figure 9: Pulse- height spectrum of dark counts measured at an overvoltage of 1.2 V and room 
temperature. The applied cuts and the high statistics of events have allowed to clearly resolved 
five peaks. A Gaussian fit was used to disentangle contributions of adjoining peaks. 



as the area under the k-th peak divided by the total number of events in the spectrum. 
Small corrections were applied to disentangle contributions of adjoining peaks in the valley 
region by means of Gaussian fitting. However, peaks are not exactly Gaussian and, conse- 
quently, differences between the experimental peak areas and those of the fitted Gaussian 
curves were added to the uncertainties, which are still dominated by statistical errors. 

In dark-count measurements, the crosstalk probability e can be estimated as the frac- 
tion of pulses with two or more triggered pixels, i.e., 1 — P(l), which was measured to be 
0.1592 db 0.0007. Accidental coincidences of non-correlated dark-counts events, however, 
would contribute to this fraction even in the absence of crosstalk. The observed dark-count 
rate, excluding afterpulses, was 1.76 MHz, and therefore, the expected average number 
of events within a time window of 3.0 ± 0.6 ns (see section 3.2) is fj, = 0.0053 ± 0.0011. 
Assuming a Poisson distribution of dark counts, the probability that a given pulse is due 
to two or more coincident dark-count events is only 0.0026. To account for this small con- 
tribution of accidental coincidences, the crosstalk probability was determined using (17), 
resulting in e = 0.1570 ± 0.0008. Note that P to t(k) = (1 - P t ot(0)) P{k) for k > 1 has to 
be used for application of equation (17), where Ptot(0) = exp(— /i). 

Theoretical probability distributions were calculated for the different crosstalk models 
explained above. For this purpose, equations (13 14) were used setting e = 0.1570 and 
assuming that P pr (m) is a Poisson distribution with mean = 0.0053. In the left-hand plot 
of figure 10 , experimental probabilities are compared with those predicted by the crosstalk 
models of [16] and [22J, based on the geometric and Borel distributions, respectively. 
The experimental data lie between these two limit situations, as expected. A geometric 
distribution (n = 1) underestimates the probability of multiple crosstalk events, while 
the Borel distribution (n — > oo) overestimates it. In the right-hand plot, the ratio of 
theoretical and experimental probabilities are displayed for these two crosstalk models 



as well as for the geometrical models described in section 2.2. Error bars represent the 



experimental uncertainties including those of the e and [i parameters used to calculate the 
theoretical probabilities. The only model that is compatible with the experimental data 
within uncertainties is the 4-nearest-neighbors model. Results of the model with n = 4 
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Figure 10: Comparison of experimental probabilities with theoretical predictions from different 
crosstalk models. Left: Experimental data lie between the two limit situations corresponding 
to both geometric and Borel distributions of crosstalk [TBI [23]. Right: Ratio of theoretical and 
experimental probabilities for different crosstalk models. Only the 4-nearest-neighbors model is 
compatible with experimental data. 



without saturation effects (section |2.1[ ), also included in the figure, deviate significantly 
from experimental data, showing that these effects are important to describe the crosstalk 
statistics. 



The excess noise factor was calculated using (21), resulting in a value of 1.1665 ± 
0.0022. Theoretical predictions of the excess noise factor for g = 0.1570 range from 
1.1570 (geometric distribution) to 1.2059 (Borel distribution). Again, the only crosstalk 
model that gives an excess noise factor compatible with the experimental one is the 4- 
nearest-neighbors model (ENF = 1.1683 is obtained numerically using the extrapolation 



of equation (12)). The comparison for the E\ and Van parameters yields to the same 
conclusion. 

Our results confirm the 4-nearest-neighbors hypothesis of crosstalk, previously used 
in Monte Carlo simulations |10H12j to describe different SiPM devices from Photonique 
and Hamamatsu. This hypothesis is further reinforced by some experiments with various 
SiPMs from Hamamatsu and SensL, performing a selective pixel illumination by means of 
a microscopic scan with a point light source 128]. The crosstalk probability was found 
to be nearly constant when illuminating inner pixels, whereas it is significantly reduced 
for pixels right at the edges of the array, since they have fewer neighbors. 

It should be emphasized that these results are not in disagreement with those of [9] 
[13], where crosstalk was found to be also induced efficiently in distant pixels in some 
non-commercial SiPM devices. In those cases, crosstalk may be better described by a 
wider neighborhood of pixels. In fact, it was shown in [22J that the experimental pulse- 
height distribution of dark counts obtained by [7j, using their own SiPM design of 1600 
pixels, fit well a Borel distribution, which is nearly equivalent to that predicted by our 
all-neighbors model. Moreover, crosstalk events with a delay of a few tens of nanoseconds 
were observed in [7j, which were attributed to a slower crosstalk process where charge 
carriers are generated in the Si bulk and migrate into the avalanche region of a perhaps 
distant pixel (see also p2]). This contribution is disregarded in our analysis, since it only 
includes simultaneous excitation of pixels within a 3 ns time window. 

4.3 Measurements under pulsed illumination 



Theoretical relationships deduced in section 2.3 were validated against experimental data 



For this purpose, a set of measurements was carried out illuminating the SiPM with the 
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Figure 11: Pulse- height spectra measured under pulsed illumination at constant light intensity 
but different bias voltages. 



O0-025 q 
> 

"50.020 ; 
O 

S 0.015 ; 

a 

&0-010 

W 

Is °- 005 
°* 0.000 




5. 



1.0 ; 

0.8 : 
0.6 : 
0.4 : 

02 : 

o.o ^A- 
o.o 



Bias voltage (V) 



A V (V) 




Figure 12: Results of the measured pulse-height spectra. Left: Distance between peaks (pro- 
portional to gain) versus bias voltage and linear extrapolation to obtain the breakdown voltage. 
Center: Average number of triggered pixels (proportional to sensitivity) versus overvoltage. Right: 
Probability of crosstalk versus overvoltage and fit of equation ( 23 ) to data. 



pulsed nitrogen laser at constant conditions (light intensity and temperature) but different 
SiPM bias voltage. Three pulse-height spectra recorded at 69.0 V, 69.4 V and 69.8 V are 



shown in figure 11 Avalanche gain is proportional to overvoltage, hence pulse amplitudes 



increase linearly at increasing bias voltage, which results in the stretching of the spectrum. 



This is illustrated in the left-hand plot of figure 12, where the distance between peaks 
(proportional to gain) is represented as a function of bias voltage. Data exhibit a pure 
linear behavior, as expected, and the extrapolation at zero gain provides the breakdown 
voltage, which is determined to be V h = 68.54 ± 0.04 stat ± 0.79 syst V at 25°C. 

Measurements were performed at very low light intensity (< 1 photon detected per 



laser shot), and a Poissonian photon statistics was assumed. Therefore, equation (16) 



allowed us to calculate the actual average number /x of detected photons from the fraction 
of laser shots in which no pulse was registered, i.e., -Ptot(O). Contribution of dark counts 
in coincidence with the laser shot (time window of 5 ns) was estimated to be less than 
1%. Results are shown in the central plot of figure 12 as a function of the overvoltage Ay, 



where data were previously normalized by the relative mean intensity of the laser beam 
measured with the photodiode during each run to compensate the laser-intensity drift. An 
approximated linear behavior is found, indicating that the SiPM sensitivity is essentially 
proportional to overvoltage (see [28J for a precise measurement of the sensitivity of the 
Hamamatsu model S10362-11-100C). 



The crosstalk probability was calculated using equation (17), where P tot (l) is deter- 



mined as described in section 4.2 Results exhibit a quadratic dependence on overvoltage 



in the experimental range (right-hand plot of figure 12), which is expected from the lin- 
ear growths of both sensitivity and gain. However, saturation should be reached at large 
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Figure 13: Comparison of experimental and theoretical results for the E\ and Var\ parameters. 
Theoretical curves are derived using the dependence of e on overvoltage given by the fit of equation 
(23) to data. 



overvoltage. Assuming that the number of crosstalk photons impinging a given neighbor 
of the primary pixel is Poisson distributed, the crosstalk probability should behave as 

s = l-exp(-KAV 2 ) , (23) 

where K is a constant of the detector. The best fit of this equation to data is also shown 
in the figure. The relatively large deviations of data at 69.0 V and 69.2 V from the fitted 
curve may be due to a lower signal-to-noise ratio introducing some bias in the calculated 
-Ptot(O) and Ptot(l) probabilities. 

Crosstalk basically contributes to shift and broaden the pulse-height spectrum as de- 



scribed in section 2.3 Experimental E\ and Var± values were obtained using (18) and 
(19), where the mean ^tot and variance Vartot were calculated for the observed proba- 
bility distribution, and E pr = Var pr = fj, is already known. The resulting E\ and Var± 



parameters are represented in figure 13 by full circles. 

Theoretical predictions for the crosstalk models are compared with these experimental 
Ei and Var\ values. Results for both the geometric and Borel distributions were calculated 
using the available analytical expressions [22], while for our proposed models they were 
numerically computed from the probability distributions of table [T] extended to k > 5 



using the geometric extrapolation of equation (12). The curves of E\ and Var\ versus 



overvoltage shown in figure 13 are derived by using the dependence of e on overvoltage 



obtained from the fit of equation (23) to data. 

All the crosstalk models yield very similar E\ and Var\ values at Ay < 1 V in 
agreement with experimental results. Note that these parameters are model independent 
up to first order in e (see ([9]) and (|10|)). However, differences at higher order become 
important at e > 0.1. In fact, predictions from models with n > 4 deviate significantly 
from experimental data at AV > 1 V. 



Unlike the measurements at dark-count conditions of section 4.2 these results do not 
allow us to clearly discriminate between crosstalk models, because the statistics is lower 
and the systematic errors are more difficult to control (e.g., photon arrival is assumed to be 
Poisson distributed with constant mean during the run). Nevertheless, our experimental 
data under pulsed-illumination are also compatible with the 4-nearest-neighbors model. 
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5 Conclusions 



The effects of optical crosstalk on the photon-counting properties of SiPMs were studied in 
detail. Novel analytical models describing the statistics of crosstalk events were developed 
and compared with experimental results. 

The proposed models take into account that each pixel of the array have a finite 
number of available neighboring pixels to excite via crosstalk. Effects of cascades due 
to crosstalk propagation from neighbor to neighbor were included in a similar way to 
the models previously reported in [16\ I22j . which can be regarded as limit situations of 
ours. In addition, geometrical considerations on the development of these cascades on the 
array, which can lead to saturation of pixels, were accounted for. As a result, analytical 
expressions for the probability distribution of the total number of triggered pixels have 
been reported assuming four different geometrical arrangements of neighbors: 4 nearest 
neighbors, 8 nearest neighbors, 8 L-connected neighbors, and all neighbors. 

Application of these analytical models to evaluate the effects of crosstalk on the photon 
statistics under pulsed-illumination conditions was also discussed. Simple expressions were 
obtained to account for the shift and broadening of the spectrum as well as for the excess 
noise factor due to crosstalk. 

A dedicated experiment, using the SiPM device S10362-11-100C from Hamamatsu, 
was carried out to validate our theoretical results. Measurements were performed at dark- 
count conditions as well as under pulsed illumination. A waveform analysis was developed 
to obtain the pulse-height spectrum with very high resolution. In particular, the influence 
of short-delay afterpulses had to be carefully taken into account. As a by-product of this 
analysis, the recovery time of pixels was estimated to be 40 ns at room temperature. 

Results of dark-count measurements were found to be consistent with the 4-nearest- 
neighbor model, in agreement with other authors that performed MC simulations of the 
response of SiPMs with crosstalk effects |10H12| . The other crosstalk models, including 
those of |16} [22j , failed to match the experimental probability distributions, which demon- 
strates that the above-mentioned geometrical considerations are important. Our results 
indicate that optical crosstalk takes place between adjacent pixels, although excitation of 
distant pixels could also be possible in a different chip configuration, as observed in [9| [T3] . 

Pulse-height spectra were measured under low-intensity pulsed illumination as a func- 
tion of the SiPM bias voltage. As expected, both the SiPM gain and sensitivity were found 
to be nearly proportional to overvoltage, and the crosstalk probability showed a quadratic 
growth. Comparison of experimental and theoretical results for pulsed illumination was 
done by evaluating the effect of crosstalk on the mean and variance of the spectrum. 
While all the considered crosstalk models can only account for experimental results at low 
overvoltage, the 4-nearest-neighbor model was found to be consistent with experimental 
data at any overvoltage, which gives a further support to results obtained at dark-count 
conditions. 
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